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基于 多 重 特征 匹配 的 点 云 配 准 算法 
李 ” 强 ， 高 保 禄 ， 窦 明亮 
(太原 理工 大 学 信息 与 计算 机 学 院 , 太原 030024) 
摘 要: el 度 低 的 问题 ， 提 出 一 种 基于 多 重 特征 匹配 的 点 云 配 
准 算法 。 首 先 采 用 改进 自 适 应 八 又 树 算 法 分 割 点 云 ， 通 过 移动 最 小 二 来 法 (MLS) 对 其 叶 节 点 进行 局 部 拟 合 后 ， 计 算 
点 的 多 重 特征 ; 然后 提出 了 基于 多 重 特征 的 点 ee 约束 的 点 对 作为 匹配 点 对 ， 进 而 求 取 旋 转 


珑 阵 和 平移 矩阵 实现 点 云 配 准 。 实 验 表 明 ， 该 算法 能 在 保持 点 云 配 准 速度 较 高 的 基础 上 ， 有 效 提升 配 准 的 准确 度 
且 准 确 度 的 提升 幅度 随 着 点 集 数量 的 增 大 呈 升 高 趋势 。 
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Point cloud registration algorithm based on multiple-feature matching 


Li Qiang, Gao Baolu, Dou Mingliang 
(College of Information & Computer, Taiyuan University of Technology, Taiyuan 030024, China) 


Abstract: To solve the problem that iterative closest point (ICP) algorithm has a single feature for searching and low 
accuracy for registration, this paper proposed a point cloud registration algorithm that based on multiple-feature matching. It 
chose the improved adaptive octree algorithm to segment the point cloud. Then calculated the multiple features of the points 
after performed moving least squares (MLS) algorithm to fit the leaf nodes. Next, this algorithm introduced the point pairs 
similarity that based on multiple features to establish the matching points. Lastly, computed the rotation matrix and 
translation matrix to achieve registration. Experiments show that this algorithm can effectively improve the accuracy of 
registration on the basis of keeping the point cloud registration speed high. And with the number of point sets increasing, the 
trend of accuracy for this method is increasing. 
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0 引言 确 匹 配点 ， 但 该 方法 在 匹配 点 选取 过 程 中 限定 的 阔 值 较 多 ; 
Elbaz 等 人 I1 引 入 深度 神经 网 络 的 方法 进行 配 准 , 但 离线 训练 
三 维 点 云 重 建 技 术 已 经 广泛 应 用 在 文物 数字 化 、 工 业 币 阶段 需要 大 量 数据 ， 更 适合 大 地 形 场景 ， 王 勇 等 人 中 通过 多 
造 建 模 等 模型 构建 领域 中 。 点 云 获 取 通 常 采 用 激光 扫描 仪 ， 分 辩 率 加 快 匹 配 速度 ， 引 入 匹配 度 概 念 改进 ICP 匹配 点 的 搜 
但 由 于 物体 的 不 可 穿 透 性 以 及 提高 准确 性 的 要 求 ， 目 标 表 面 。” 索 ， 但 是 该 方法 匹配 点 的 特征 单一 ， 在 点 集 较 大 时 会 产生 错 
完整 信息 往往 需要 经 过 多 和 角度、 多 次 扫描 获得 ， 导 致 获取 的 。” 误 匹 配点 对 ， 准 确 度 一 般 。 
点 云 数 据 不 在 同一 坐标 系 下 ， 因 此 需要 对 不 同 角 度 、 不 同 批 本 文 提 出 了 一 种 基于 多 重 特征 匹配 的 点 云 配 准 算法 。 首 
次 获取 的 点 云 进 行 配 准 。 点 云 配 准 的 实质 是 求解 待 配 准点 云 先 ， 通 过 改进 自 适 应 八 又 树 算法 将 点 云 数 据 组 织 在 空间 单元 
到 目标 点 云 的 刚性 转换 关系 。 中 ， 对 每 个 单元 采用 移动 最 小 二 乘法 (moving least squares， 
迄今 为 止 ， 针 对 点 云 配 准 问题 已 提出 了 很 多 方法 ，Besl ”MLS) 拟 合 局 部 曲面 后 计算 出 点 的 多 重 特征 ， 然 后 利用 多 重 
等 四 提出 的 最 近 点 友 代 算法 ， 即 ICP(iterative closest point) ， 特征 代替 文献 [7] 中 单一 的 曲率 特征 ,并 采用 基于 多 重 特征 的 
是 应 用 最 广 的 一 种 。 该 算法 取 目 标点 云 与 待 配 准 点 云 间 欧 氏 点 对 相似 度 代 蔡 其 中 的 匹配 度 建 立 匹 配点 对 ， 最 后 用 四 元 数 
距离 最 小 的 点 来 构建 匹配 点 对 ， 求 得 匹配 点 对 的 旋转 和 平移 ”法 完成 点 云 配 准 。 


变换 和 矩阵， 变换 后 重新 建立 匹配 点 对 进行 迭代 计算 ， 达 到 收 
敛 约 束 条 件 后 完成 配 准 。 该 算法 构建 匹配 点 对 的 过 程 计算 量 1 改进 自 适 应 八 又 树 算法 

非常 大 , 且 特 征 不 明显 区 域 元 余 的 点 大 幅 降低 了 其 收敛 效率 。 ”1.1 自 适应 八 又 树 算 法 

针对 其 缺陷 ， 国 内 外 学 者 提出 了 众多 改进 ICP 算法 。Sharp 点 云 分 割 中 八 又 树 的 划分 策略 直接 影响 着 划分 结果 所 占 
等 人 馈 以 点 到 邻 域 重心 的 距离 作为 特征 进行 点 云 配 准 ， 以 重 ”的 存储 空间 和 构建 时 间 。 自 适应 八 又 树 通过 设置 阔 值 ， 可 以 


出 


心 描述 单个 点 在 点 云 中 的 相对 位 置 ， 但 由 于 特征 单一 ， 准 确 ”将 每 个 节点 中 点 的 数量 控制 在 合理 范围 内 ， 同 时 减少 分 割 次 
度 提 升 较 小 ， 韦 盛 斌 等 人 外 中 以 点 面 距 离 构 建 匹 配点 对 ,为 数 ， 是 点 云 最 常用 的 分 割 算 法 之 一 ， 但 是 其 闵 值 的 设 定 不 够 
点 云 特征 的 选择 提出 了 新 思路 ， 但 由 于 曲面 的 拟 合 较 差 影响 。 灵活 ， 如 果 闵 值 太 大 , 无 法 体现 分 割 的 效果 ; 如 果 闪 值 太 小 ， 
了 准确 度 ， 杨 小 青 等 人 铝 采 用 主 曲率 约束 和 点 云 法 向 量 夹 角 每 个 节点 内 的 点 太 少 ， 根 据 最 小 二 乘法 的 原理 一 一 通过 最 小 
初步 选取 点 集 ， 然 后 采用 点 间距 离 和 高 斯 曲率 进一步 获取 精 “化 误差 的 平方 和 寻找 数据 的 最 佳 函数 匹配 , 由 于 选取 的 点 少 ， 
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录用 定稿 李 强 ， 等 : 基于 多 重 特征 匹配 的 点 云 配 准 算法 第 37 卷 第 2 期 
放大 了 边缘 点 、 离 群 点 以 及 噪声 的 误差 对 于 曲面 生成 的 影响 ， 直到 叶 节 点 中 的 数目 小 于 Daax， 然 后 在 叶 节 点 中 ， 采 用 PCA 
降低 了 拟 合 曲面 准确 度 ， 为 了 便于 理解 ， 图 1(a) (b) 以 二 维 算法 计算 点 法 向 量 信息 ， 进 而 求 得 法 向 量 标准 偏差 ag， 给 定 
线 拟 合 说 明了 离 群 点 对 不 同 数据 集 的 影响 ， 对 于 图 中 (a) 和  ” 法 向 量 标准 偏差 的 阔 值 &, 若 o>&, 则 说 明 曲 面 在 该 节点 曲率 
(b)， 在 采用 移动 最 小 二 乘法 拟 合 时 ， 若 都 增加 离 群 点 A， 其 变化 比较 大 ， 须 继续 细 分 ， 直 到 叶 节 点 中 的 数目 小 于 Dain; 
对 (a) 中 拟 合 曲线 的 影响 将 会 明显 大 于 (b)。 同 理 ， 对 于 复杂 的 否则 不 再 细 分 ， 最 后 采用 MLS 算法 对 叶 节 点 中 的 数据 进行 
点 云 模型 中 ， 被 分 割 物体 表面 常常 会 有 连续 的 平顺 区 域 〈 如 曲面 拟 合 ， 求 得 其 局 部 拟 合 曲面 z=flx,y)。 
图 1(c) 中 和 拖 形 框 区 域 ) ， 若 靖 值 太 小 ， 分 割 后 每 个 块 中 的 点 给 定点 集 G 的 法 向 量 的 标准 偏差 c 的 计算 公式 为 式 (1)， 
数量 太 少 ， 反 而 会 降低 后 续 MLS 拟 合 精确 度 。 反映 了 点 集 G 内 曲面 曲率 变化 的 剧烈 程度 四。 
1.2 改进 自 适 应 八 又 树 FI Rr pr 2 

本 文 提 出 了 改进 自 适 应 八 又 树 ， 首 先 采用 自 适 应 八 又 树 .- 恬 (1) 
算法 ， 设 置 较 大 的 阔 值 对 点 云 数据 粗 分 制 ， 在 此 基础 上 提出 Ne 
了 增加 法 向 量 的 标准 偏差 判断 每 个 八 又 树 节点 曲面 的 变化 剧 经 过 分 析 ， 改 进 自 适应 八 又 树 ， 一 方面 减少 了 不 必要 的 
烈 程度 ， 如 果 大 于 给 定 闵 值 就 继续 细 分 ， 直 到 满足 叶 节 点 的 ” 分割 从 而 减少 了 后 续 拟 合 次 数 ， 另 一 方面 由 于 采用 了 法 向 量 
最 小 冰 值 ， 最 后 对 分 割 后 的 每 个 时 节点 进行 局 部 拟 合 ， 考 虑 。 标准 偏差 判断 ， 有 效 降 低 了 每 个 八 叉 树 叶 节 点 内 点 云 拟 合 
到 本 文 的 目的 是 对 点 云 数据 进行 配 准 ， 因 而 可 以 避 开 因此 带 机 的 复杂 度 ， 所 以 有 助 于 后 续 点 云 处 理 中 局 部 曲面 的 快速 精 
来 的 曲面 拼接 问题 由 。 确 拟 合 。 

分 割 过 程 如 下 : 设 定 叶 节点 数量 的 最 大 最 小 值 分 别 为 算法 的 形式 化 说 明 如 图 2 所 示 。 

分 割 阶段 ， 采 用 递归 的 方式 进行 循环 分 割 ， 
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(a) 少 量 点 拟 合 
图 1 
Fig. 1 


New Octree: 
New memory; 


Create list of pointers; 


New OctreeRoot; 


Subdivide_Node: 
Subdivide 

for 遍历 8 个 子 节点 

| 芷 节点 内 非 空 


/创建 新 八 又 树 

/申请 足够 大 的 内 存 
/创建 指针 列表 

/创建 八 叉 树 根 节点 


/改进 的 自 适 应 八 又 树 


/将 传 入 节点 分 成 8 个 


(b) 大 量 点 拟 合 
离 群 点 对 MLS 的 影响 


Influence of outliers on MLS 


ff 节点 内 点 的 数量 >D,。、 
Subdivide_Node ”// 递 归 细 分 节点 


在 点 集 区 域 Q 上 建立 拟 合 函数 00， 表 示 为 式 (2) 。 
f= a= Wa O) 
qg(0)=[L x, y,%,, 7 (3) 
b=(D#) (4) 
其 中 : 所 求 系数 a(x)=[az(X),a2(X?),A,am(X)]7 为 x 的 函数 ; 基 函 


数 400=[9(00,4g200,49gnC0] 为 上 阶 完备 多 项 式 ， 式 〈3) 为 


| 
else 下 ( 节 点 内 点 的 数量 二 D，，、 y 离 
| && 节点 内 点 的 数量 >D，) 其 二 次 基 ， 取 m=6; 式 (4) 为 向 量 x=[xi,x2,4,xn] 的 加 权 离 
| | | 计算 节点 内 的 法 向 标准 偏差 5 散 LL 范式 。 
| | | fo>é /标准 偏差 小 于 阔 值 E 
| | | | Subdivide_ Node /递归 细 分 节点 J=》> wz-z)LACO-y> 下 = 
| | | else 六 (5) 
| I > we nat] 
Continue 2 二 a yw 
式 (5) 中 4 为 区 域 Q 内 点 的 数目 ，f) 为 拟 合 函数 ， 
| Delete /删除 当前 节点 处 的 权 函 数 为 w(x-xj) ，a(x) 为 待定 系数 。 求 导 可 得 
end; /结束 循环 Oo 
=A(Wa(W) -By=0 6 
图 2 算法 的 形式 化 说 明 ee 
Fig.2 Formal description of algorithms Q(x)=A CD)BCDy (7) 
日 . CH 
2 ”点 云 处 理 , 
4CD= w(x—x,)g(x)g (x)) 8 
2.1 MLS 拟 合 > (8) 
本 文选 用 MLS 法 对 分 割 后 的 叶 节 点 数据 进行 曲面 拟 合 ， B(x) =[w(x-X)g(m), w(x-%)g(%,), (9) 
求 得 其 局 部 拟 合 曲面 z=ftx,y)， 算 法 过 程 如 下 : A, w(x-%,)g(%)] 
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此 时 即使 g(x) 为 多 项 式 , 式 (11) 中 的 fw) 也 不 再 是 多 项 
式 。 如 果 基 函数 g EC' 则 权 函 数 w EGC, 则 拟 合 函数 fEC™"("5) 
2.2 多 重 特征 计算 
在 对 曲面 的 局 部 几何 特征 进行 量化 比较 时 ， 由 于 待 配 准 
据 是 同一 个 物体 在 不 同 坐 标 系 下 的 两 组 点 云 ， 因 此 所 选 的 
特征 就 必须 具有 缩放 、 旋 转 以 及 平移 不 变性 的 特征 ， 经 过 得 
选 ， 本 文 使 用 高 斯 曲率 及、 平均 曲率 瓦 、 主 曲率 kl 和 kz、 与 
重心 的 距离 L、k 近邻 点 法 向 量 夹 角 的 平均 值 M 来 描述 点 云 
中 某 一 点 所 在 曲面 的 局 部 特征 。 
其 中 ， 高 斯 曲率 K 可 以 衡量 曲面 在 该 点 处 总 的 弯 

_LN-M’ 

EG-F’ 
面 在 该 点 的 平均 弯 
EN-2FM +GL 
2(EG-F’) 


以 用 来 对 


程度 为 
Q4) 


平均 曲率 H 表示 程度 为 


hs (15) 


面 在 该 点 处 指定 方向 


两 个 主 曲 率 ki 和 ko 可 
的 弯曲 程度 进行 度量 为 


kh=H -NH’-K 


(16) 


k=H + NH’-K 
与 重心 的 距离 L 描述 了 该 点 在 点 云 中 的 相对 位 置 。 
(2 C2] + -> | 


点 集 P 中 的 任意 其 近邻 点 法 向 量 夹 
值 Mi 反映 了 点 云 局 部 曲面 的 尖锐 程度 。 


1 元 :元 
M, = 天 Zarccos 一 | 
疝 


元 | 
综 上 所 述 ， 对 于 点 集 P 中 的 任意 一 点 p:， 根 据 式 (10) 


得 到 了 其 所 在 叶 节点 的 局 部 拟 合 曲面 z=ftx,y)， 再 由 式 (14) 
~《〈19) 取得 了 其 局 部 多 重 特征 。 


3 ”基于 多 重 特征 匹配 的 点 云 配 准 


最 近 点 迭代 算法 (ICP) 采用 最 近 点 匹配 策略 ， 缺 点 是 
易 产生 大 量 错误 的 匹配 点 对 ， 特 别 是 待 配 准 点 云 与 目标 点 


(17) 


(18) 
的 平均 


一 点 Di 


(19) 


经 典 ICP 算法 ， 是 对 待 配 准 上 
目标 点 云 (GP) 中 每 个 点 的 距离 ， 选 取 
最 近 点 作为 匹配 点 ， 这 样 就 建立 了 从 BP 到 GP 
式 (20) 计算 R 和 T， 但 因 
所 以 需要 重复 运行 


计算 出 其 


(GP) 中 的 
应 关系 ， 然 后 采用 


的 一 一 对 
样 建立 的 


函数 的 最 优 解 ， 利 月 
旋转 矩阵 R 和 平移 矩 
距离 平方 和 作为 目标 函 


f(H)=f (RT)=2, 
点 集 为 Pi， 


示 的 点 间 


本 文 
BP 和 GP 
匹配 点 对 


首先 采 | 
样 ， 该 方法 先 采 档 
点 提升 配 准 精度 ， 最 后 小 
式 (22) 得 到 的 向 旨 
最 大 分 辨 率 为 n， 帆 
m) 级 提取 点 的 采 档 


样 根据 


合 
云 重 且 度 较 小 的 情况 下 ,文献 [7] 的 改进 算法 采用 的 曲率 匹配 
度 仍 存在 特征 单一 、 准 确 度 较 差 的 问题 。 本 文 基于 多 重 特征 
的 概念 ， 在 增加 约束 特征 的 同时 采用 了 更 
加 适合 多 重 特征 匹配 的 向 量 相 似 度 作 为 判断 条 件 。 有 效 提 高 
准确 性 。 
点 云 的 配 准 原 理 就 是 求 待 配 准 云 向 目标 点 云 转换 的 旋转 
和 矩阵 R 和 平移 矩阵 的 过 程 ， 如 图 3 所 示 ， 假 设 左上 方 GP 
(gray points) 是 目标 点 云 ， 右 下 方 BP (black points) 是 待 
配 准 点 云 , 该 算法 就 是 计算 BP 到 GP 的 平移 和 旋转 矩阵 ,使 
BP 和 GP 尽量 习 


量 重 又 。 


Pp: 待 配 准 


其 中 ; countm》 


与 


对 应 关系 是 
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中 ,通过 


的 准 


小 用 


为 向 零 取 整 。 


然后 ， 
多 重 特征 


对 


主 曲率 ki、 
pis 为 其 
取 点 pi 在 


特征 向 量 (qj/，4q 


弦 相 似 度 ， 


912， 


公式 为 


了 


图 3 配 疹 
Fig.3 Registration 
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y” =[y,,y2,A,y,] (10) 其 旋转 过 程 如 式 (20〉 所 示 。 
将 式 (7) 代入 式 (2) 即 可 得 出 移动 最 小 二 乘 拟 合 函数 ， 2 2 
f= = 00y GD 时 全 加 
其 中 : k 为 基 函 数 的 阶 数 ，Cx9 是 形 函 数 根据 上 述 配 准 原理 ， 若 求 取 旋转 矩阵 R 和 平移 矩阵 T 
Ox) =[@to2,A,@] (12) 就 需要 获得 点 云 数 据 BP 和 GP 中 点 的 一 一 对 应 关系 ,但 是 在 
=g" (x)A(x) Bx) 点 云 数 据 获 取 过 程 中 ， 三 维 激光 扫描 仪 是 通过 激光 测 距 的 原 
如 果 k=0， 则 基 函 数 g(x)={1}， 这 时 形 函 数 为 式 (13) 理 记录 被 测 物体 表面 大 量 的 密集 点 的 三 维 坐 标 ， 具 有 一 定 随 
的 Shepard 函数 : Wa a 角度 等 的 影响 ， 获 得 的 点 云 数 据 并 没有 
Opard (x) = W(X 一 Xi ) 绝 天 
= 


-0.5 


E 示 


意图 


点 


直到 均 方 差 误差 小 于 给 定 阔 


个 假设 ， 


schematic diagram 


i 云 上 的 (BP) 每 个 点 ， 先 


值 。 


目标 点 云 


计算 过 程 ， 


故 ICP 算法 实质 是 求解 目标 


巨 阵 TI 


]5 


日 最 小 二 乘 原则 计算 得 到 使 目标 函数 最 小 


目标 函数 常 选取 如 
数 021， 


式 (21) 所 


点 对 相 
确 度 ， 减 少 了 错误 
多 分 辨 率 关 键 点 采样 法 中 对 点 云 数 据 进 行 采 
# 少 数 点 进行 迅速 配 准 ， 
局 增加 采样 ， 
量 夹 角 将 点 云 中 的 点 分 为 m 级 ， 设 


Count, . 


以 度 最 


下 


lo -RP-7Th 


(21) 


标点 集 为 Qi, N 是 匹配 点 对 总 数 。 
基于 多 重 特征 提出 了 点 对 相似 度 概念 ， 
高 的 建立 匹配 点 对 ， 提 


在 上 述 点 云 
高 了 


配 。 


然后 大 幅 增 加 采样 
全 点 云 整体 收 剑 。 采 


1 分 辩 率 为 1 (1 三 1<n) 时 , 第 s (1 入 5 二 
比例 为 
1 


—mts 


count, l+e 


第 m 级 总 占 点 数 ，count; 为 第 s 级 总 点 数 , fix 


于 P 中 的 任 


9Gj13， 


-fix(m 


光 (22) 


/2)H) 


k2, pi3、 Pi4 分 别 为 
到 重心 的 距离 L, pi 为 该 点 法 向 量 夹 角 平均 值 M。 求 
标点 集 Q 对 应 的 k 近邻 点 gj， 


gi) 


Gj4y 


意 点 pi， 基 于 点 的 多 重 特征 建 并 其 
向 量 (pijy，piz，pi3，pi4，pis，pi6)， 其 中 pii、 


p22 为 其 


高 斯 曲率 、 平 均 | 


率 H， 


qdi5, 


辣 


同样 建立 qj 的 多 重 
其 相似 度 采用 向 量 余 
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sf 4 ”实验 结果 及 分 析 
LA (23) PN 3 、 
4 本 文 使 用 Matlab2016a 对 算法 进行 编程 ， 操 作 系统 为 
加 Windows7、CPU 为 Intel Core i7、 内 存 8 GB 、 显 卡 为 NVDIA 
分 别 计算 pi 与 其 k 近邻 点 gj; 的 相似 度 W(pi:，gq)， 以 W GeForce GT650M 并 配 有 256 GB 固态 硬盘 。 实验 数据 来 自 斯 
值 最 大 的 点 作为 pi 的 匹配 点 。 坦 福 大 学 3D 点 云 数据 库 中 的 Dragon 和 Bunny 点 云 数 据 。 
最 后 采用 四 元 数 法 I 计算 匹配 点 对 间 的 旋转 矩阵 R 和 ”4.1 不 同 阐 值 八 又 树 分 割 后 的 曲面 拟 合 速率 比较 
平移 矩阵 T。 为 了 选取 改进 和 本 文 进行 
接着 对 算法 的 空间 和 时 间 复 杂 度 进行 分 析 , 在 相似 度 W 了 大 量 重复 实验 ， 实 验 选 取 模 型 为 Bunny 模型 ， 共 33279 个 
的 计算 中 ， 其 增长 最 快 的 项 是 二 次 方 ， 因 而 相似 度 计 算 的 时  ” 点, 采用 / Te 列 如 Dmax 首次 取 总 点 数 1/8 的 
间 复 杂 度 为 O(n*)， 而 对 于 P 中 的 每 个 点 pi 与 都 要 与 其 k 近 ”大约 值 4000， 然 后 采用 减 半 法 逐次 递减 ，Dmin 取 值 为 Dmax 
邻 点 人 计算 一 次 相似 度 ，P 中 的 mm 个 点 就 需要 计算 mik 次 ， 的 /18，&=0.0001 归 ， 实 验 结果 如 表 1 所 示 。 为 了 对 下 一 步 的 
寻 此 算法 的 时 间 复 杂 度 为 Oo3); 算法 计算 过 程 中 每 个 点 的 需 。 点 云 处 理 提供 参考 ， 本 文 耗 时 选取 的 是 基于 改进 自 适应 八 又 
要 存储 的 数据 有 多 重 特征 和 具有 最 大 相似 度 的 点 ， 因 此 算法 树 的 局 部 MLS 拟 合 总 时 间 。 
所 需 的 存储 单元 是 随 着 点 的 数量 成 线性 增长 ， 即 算法 的 空间 表 1 列举 了 不 同 的 Dmax，Dmin 取 值 对 应 的 耗 时 (单位 : 
复杂 度 为 O(n)。 秒 ) ， 从 中 可 以 看 出 在 所 选 数据 范围 的 结果 中 ，Demax=1000， 
é 开始 六 Dmin=125 时 总 耗 时 最 小 ， 并 且 无 论 是 变换 Daax 还 是 Dmin 算 
本 y 一 法 总 耗 时 朝 各 个 方向 呈 增 大 趋势 。 经 分 析 是 Daax 增 大 以 后 虽 
计算 目标 点 云 p; 的 多 重 特征 然 减少 了 分 割 次 数 s， 但 是 增 大 了 局 部 拟 合 的 难度 ，Dnmin 减 
T 小 以 后 虽然 降低 了 局 部 拟 合 难度 ， 但 是 是 拟 合 次 数 增多 ， 增 
PER 加 了 系统 资源 的 调度 消耗 ， 因 此 也 增 大 了 总 耗 时 。 
i be da 在 此 基础 上 经 过 多 次 实验 ， 最 终 选 择 了 总 耗 时 更 小 的 参 
数 ， 令 Dax=1000，Dmin=100， 其 总 耗 时 为 10.312s。 
根据 分 辩 率 对 P 中 数据 进行 采样 表 1 算法 拟 合 过 程 比较 1 
y Table 1 Comparisonsl of fitting processes of algorithms 
求 取 p 在 目标 点 云 Q 中 对 应 的 匹配 点 Do i 
4000 2000 1000 500 
! 500 18.959 17.527 16.313 无 
通过 四 元 数 法 求 取 旋 转 垂 阵 R 和 平移 矩阵 T 250 18.514 15.868 13.973 15.969 
下 y 耕 125 14.753 13.031 11.711 13.632 
通过 R 和 T 对 点 集 P 变 换 后 得 到 P?” 62 16.319 14.918 13.716 15.127 
v 4.2 曲面 拟 合算 法 比较 
款 冰 数值 是 否 满足 条 件 为 了 验证 基于 改进 自 适应 八 又 树 的 局 部 MLS 拟 合 算 
法 的 速度 和 准确 度 ， 分 别 与 移动 最 小 二 乘法 拟 合 曲面 算法 、 
下 一 文献 [14，15] 中 的 算法 进行 对 比 。 选 择 相 同 的 点 云 数 据 ( 均 
- _ 分辩 率 是 Xn 采用 Bunny 模型 ) ， 通 过 多 次 重复 实验 ， 选 取 Dmax=1000， 
En Dmin=100，&=0.0001， 表 2 通过 分 块 数量 、 分 割 耗 时 、 拟 合 
_ 耗 时 、 算 法 总 耗 时 对 算法 的 速度 进行 了 分 析 ， 残 差 值 取 曲 国 
> = 2 残 差 n06' 7 的 绝对 值 均值 ， 值 越 低 越 好 。 
图 4 算法 流程 图 由 表 2 可 以 看 出 ， 直 接 对 未 处 理 点 云 采 用 移动 最 小 二 乘 
Fig.4 Algorithm flow 拟 合 ， 不 仅 耗 时 最 长 且 残 差 值 也 最 高 ， 文 献 [14] 采 用 了 基于 
基于 多 重 特征 匹配 的 点 云 配 准 算法 流程 图 见 图 4， 其 详 自 适应 八 又 树 的 MLS 拟 合 ， 由 于 分 块 数量 多 ， 拟 合 时 间 比 
细 步 又 如 下 : 较 长 ， 文 献 [15] 对 点 云 分 割 后 选取 的 拟 合 方法 是 紧 文 撑 径 向 
a) 设 待 配 准点 集 为 P， 目 标点 集 为 Q。 计 算 P 中 所 有 点 ， 基 函 数 ， 该 方法 拟 合 曲面 光滑 ， 残 差 值 最 低 ， 但 是 也 因为 待 
Pi 的 多 重 特征 ; 求解 方程 组 维 数 较 大 ， 运 算 量 大 ， 增 大 了 时 间 开 销 。 本 文 算 
b) 根 据 法 向 量 夹 角 平 均值 Mi 将 P 中 的 点 分 为 m 级 并 设 ”法 由 于 增加 了 法 向 量 标准 差 判断 叶 节点 的 特征 复杂 性 ， 有 效 
置 最 大 分 辩 率 n; 减少 了 八 又 树 的 分 割 ， 对 比 文献 [14]， 更 低 的 残 差 说 明了 更 
c) 采 用 多 分 辩 率 关键 点 采样 法 对 处 理 点 云 ， 获 得 分 辨 率 ”合理 的 分 割 ， 对 比 文献 [15]， 在 拟 合 准确 性 基本 相同 的 情况 
t 和 采样 比例 Rs; 下 ， 大 幅 降 低 了 算法 总 耗 时 。 
d) 由 式 (23) 求 取 采 样 点 的 匹配 点 ; 图 5 是 基于 四 种 方法 绘制 的 残 差 图 , 为 了 提高 展示 效果 ， 
e) 对 步骤 d) 获 得 的 匹配 点 对 使 用 四 元 数 法 计算 ， 获 取 旋 图 中 随机 选取 1400 个 点 ,结合 点 的 分 布 和 纵 坐 标的 范围 可 以 
转 矩 阵 R 和 平移 矩阵 T; 看 出 ， 文 献 [15] 的 残 差分 布 最 接近 x 轴 ， 且 纵 坐 标 范围 最 小 ， 
f) 使 用 步骤 e) 获 得 的 R 和 TT 对 了 利用 式 (20) 进 行 变换 ， 说 明 拟 合 的 准确 性 最 好 ， 本 文 算法 由 于 加 入 了 法 向 量 偏 差分 
得 到 新 的 点 云集 P; 析 ， 因 此 和 文献 [15] 几 乎 相同 ， 文 献 [14] 残 差分 布 范 围 较 广 ， 
是 重复 步骤 d)~f)) 直 到 满足 目标 函数 ; 说 明 其 拟 合 的 曲面 变化 较 大 ， 而 移动 最 小 二 乘法 的 残 差 几乎 
h) 若 分 状 率 不 为 n"， 则 进入 下 一 分 状 率 ， 返 回 步 又 c)。 随机 分 布 ， 且 范围 较 大 ， 表 明 点 云 近 似 均匀 的 分 布 在 曲面 两 
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侧 ， 并 未 拟 合成 有 效 的 表面 模型 。 经 典 ICP 算法 、 文 献 [7] 的 改进 ICP 算法 实验 结果 进行 对 比 。 
表 2 算法 拟 合 过 程 比较 2 由 于 两 组 点 云 无 一 一 对 应 关系 ， 业 内 对 评价 配 准 误差 的 方法 
Table 2 ” ”Comparisons2 of fitting processes of algorithms 还 未 形成 普遍 认可 的 方法 ， 本 文采 取 文 献 [13] 中 的 方法 来 评 
分 块 数量 /个 “分割 耗 时 /8 拟 合 耗 时 /s 算法 总 耗 时 4 残 差 值 介 配 准 误差 。Bunny 模型 的 配 准 实验 结果 如 图 6 所 示 ， 其 中 
移动 最 小 二 图 6(a) 为 配 准 前 的 Bunny 模型 ,来 自 两 个 视角 (0 度 和 45 度 ) 
乘 拟 合 全 YN 并 存在 平移 现象 ， 点 数 分 别 为 36580 和 33279， 图 6 (b) 为 
文献 [14] 337 5.667 26.614 28.281 0.02180 经 典 ICP 算法 ， 可 以 看 见 在 耳 人 条 、 脚 趾 部 分 由 于 特征 复杂 ， 
文献 [15] 337 5.596 37.156 38.752 0.00923 偏差 较 大 ; 图 6 (c) 为 文献 [7] 算 法 的 配 准 结果 ， 在 耳 杀 部 分 
本 文 算法 89 0.727 9.585 10.312 ”0.01031 较 经 典 ICP 有 明显 改进 ,但 是 特征 复杂 区 域 点 仍 有 明显 偏差 ， 
其 是 边缘 较为 杂乱 ; 图 6(d) 为 本 文 所 提出 的 算法 配 准 结果 ， 
0.4 经 过 实验 点 对 相似 度 取 99.9%， 可 以 看 出 边缘 清晰 ， 配 准 效 
0 : 果 强 于 上 述 两 种 算法 。 
0.2 全 i 
人 > ».: 


人 ) 文 献 [14] 算 法 


图 6 Bunny 模型 的 配 准 实验 结果 


表 3 三 种 算法 配 准 结果 对 比 


ICP 算法 


(c) 文 献 [7] 算 法 (dd) 本文 算法 


Fig.6 Registration experimental results of Bunny model 


Table 3 Comparison of registration results of three algorithms 


配 准 方法 配 准 耗 时 /s 配 准 误差 

经 典 ICP 算法 120.1214 0.0048 

(c) 文 献 [15] 算 法 文献 [7] 13.1456 0.0025 

0.06 本 文 算法 11.8748 0.0014 
. . ee 表 3 为 三 种 算法 在 配 准 速度 和 配 准 误差 上 的 横向 对 比 ， 
we 经 典 ICP 算法 耗 时 最 多 ， 主 要 原因 是 匹配 点 的 搜索 和 进 代 过 
A 程 较 慢 ; 文献 [7] 算 法 采用 了 多 分 辨 率 和 kd-tree 加 速 ， 并 通过 
> 二 曲率 特征 选择 匹配 点 ， 在 准确 度 和 速度 上 都 有 一 定 提升 ， 而 
本 文 提出 的 基于 多 重 特征 匹配 的 点 云 配 准 算法 ， 虽 然 在 求 取 
Se 点 多 重 特征 的 过 程 中 耗 时 较 多 ， 但 是 通过 改进 的 自 适应 八 叉 
“3 We 树 大 幅 减少 了 分 割 块 和 局 部 拟 合 的 次 数 ， 反 而 配 准 耗 时 略 有 
-0.06 缩短 。 在 配 准 误差 上 ， 一 方面 合理 的 分 割 提高 了 多 重 特征 的 
(qd) 本文 算法 精确 性 ， 另 一 方面 基于 多 重 特征 的 点 对 相似 度 使 得 点 匹配 更 

图 5 四 种 方法 残 差 图 加 精确 ， 从 而 配 准 误差 有 了 明显 的 下 降 。 
Fig.5 Residual diagrams of four method 为 了 对 本 算法 的 有 效 性 做 进一步 验证 ， 使 用 更 加 复杂 的 
4.3 点 云 配 准 算法 实验 点 云 数 据 Dragon 实验 ， 图 7(a) 为 Dragon 配 准 前 的 点 云图 ， 
局 部 曲面 拟 合 是 为 了 获得 点 的 多 重 特征 ， 本 文 算法 最 终 。 ”点 数 分 别 为 435545 和 418387。 

目的 是 点 云 配 准 。 为 了 验证 本 文 所 提 算 法 的 有 效 性 ， 分 别 和 由 于 点 多 图 小 ， 图 7 中 展示 效果 不 理想 ， 因 此 对 其 局 部 


放大 进行 说 明 ， 


李 强 ， 等 : 


图 8 是 龙头 部 分 的 放大 图 。 从 图 8 中 可 以 看 


ICP 配 准 ， 龙 f 


分 点 集 杂 乱 无 音 


; 文献 [7] 算 法 配 准 效果 如 图 


所 改进 ， 但 是 尖锐 部 分 的 离散 率 
用 了 曲率 作为 匹配 依据 ， 导 致 在 点 密度 大 、 特 征 
误 率 仍然 较 高 ， 而 本 文 算 法 在 相似 度 取 值 
如 图 9(d)， 边 缘 清 晰 ， 点 和 


中 匹配 点 的 错 


99.99% 时 配 准 效果 


9(c)， 
万 然 较 大 ， 这 是 由 于 


图 8 中 龙 角 尖锐 部 分 的 放大 图 ， 从 图 9(b) 可 以 看 出 采用 经 典 
边缘 点 离散 率 较 大 ， 边 缘 模 糊 不 清 ， 


图 9 是 
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算法 增幅 略 大 , 这 是 由 于 本 文 算法 虽然 在 分 割 阶段 耗 时 较 少 ， 


但 是 由 于 配 准 过 程 中 点 的 特征 的 计算 较 多 ， 匹 配 阶 段 基于 多 


重 特征 的 点 对 相似 度 计算 也 增加 了 算法 耗 时 。 而 在 配 准 误差 


中 心 部 


其 边缘 有 


口 7 
从 个 水 


低 了 匹配 错误 率 。 


(a) 原 始 数 据 (b) 经 


较为 规整 ， 配 
准 效果 最 好 ， 这 是 因为 采用 了 多 重 特征 匹配 ， 多 特 和 


复杂 的 区 域 


FE 约 束 降 


ICP 算法 


(c) 文 献 [7] 算 法 


(a) 原 始 数 据 (b) 经 


图 7 Dragon 模型 的 配 准 实验 结 呈 


R 


(d) 本 文 算法 


Fig.7 Registration experimental results of Dragon model 


(c) 文 献 [7] 算 法 (d) 本 文 算法 


图 8 ”龙头 配 准 放大 


Fig.8 Registration enlargement map of crown 


表 4 三 种 算法 配 准 结果 对 比 


Table 4 Comparison of registration results of three algorithms 


误差 增幅 ( 倍 ) 


配 准 方法 配 准 耗 时 /s 增幅 ( 倍 ) 配 准 ; 
经 典 ICP 算法 1740.4877 14.49 0.1076 22.42 
文献 [7] 121.6585 9.25 0.0351 14.04 
本 文 算法 117.4356 9.89 0.0173 12.36 
表 4 说 明了 随 着 点 云 数量 的 大 幅 增 加 ， 三 种 方法 耗 时 也 
增幅 明显 ， 其 中 传统 ICP 算法 增幅 最 大 ， 本 文 算法 较 文献 [7] 


本 文 算法 增幅 最 小 ， 进 一 步 说 明了 本 文 的 算法 在 点 云 配 
Fh 具有 较 高 的 准确 度 。 


配 准 耗 时 长 、 准 确 度 低 的 问题 ， 提 出 了 一 种 基于 多 重 特征 匹 
配 的 点 云 配 准 算法 。 首 先 提 出 了 改进 自 适 应 八 又 树 算法 分 割 
点 云 ， 以 MLS 作为 曲面 拟 合算 法 对 分 割 后 的 点 云 数据 生成 


(a) 原 始 数据 (b) 经 典 ICP 算法 
Y : 2 2 x -> 六 $ 
(c) 文 献 [7] 算 法 (d) 本 文 算法 


图 9 龙 角 配 准 放大 图 


Fig.9 Registration enlargement map of Longjiao 
结束 语 
本 文 针对 点 云 配 准 过 程 中 匹配 点 对 搜索 规则 单一 导致 的 


局 部 拟 合 曲面 ， 并 基于 此 计算 点 的 多 重 特征 ， 然 后 提出 了 基 


于 多 重 特征 的 点 对 相似 度 概念 ， 利 用 点 对 相似 度 得 到 更 加 精 
确 的 匹配 点 对 ， 实 现 点 云 配 准 。 


扩大 了 曲面 拟 合 的 数据 量 , 降低 了 噪声 点 、 离 群 点 等 的 影响 ， 


实验 表明 ， 本 文 方法 有 效 减 少 点 云 分 割 块 数 的 同时 合理 


而 提高 了 点 多 重 特征 的 准确 性 ; 配 准 阶段 基于 多 重 特征 的 
对 相似 度 大 大 减少 了 错误 匹配 点 ， 提 高 了 配 准 准 确 度 。 
1 


晶 是 ， 本 文 方法 仍 存在 诸多 不 足 之 处 ， 一 是 在 分 割 点 云 


时 闹 值 需 根据 点 云 的 大 小 手动 调 参 ， 较 为 繁琐 ， 二 是 采用 多 


重 特征 ， 增加 了 每 个 点 需要 存储 的 信息 , 增 大 了 空间 复杂 度 ， 
这 些 


些 是 本 文 在 接 下 来 的 研究 中 将 要 进一步 完善 的 方向 。 
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